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ABSTRACT 



Context. The winds from a non-accreting pulsar and a massive star in a binary system collide forming a bow-shaped shock 
structure. The Coriolis force induced by orbital motion deflects the shocked flows, strongly affecting their dynamics. 
Aims. We study the evolution of the shocked stellar and pulsar winds on scales in which the orbital motion is important. 
Potential sites of non-thermal activity are investigated. 

Methods. Relativistic hydrodynamical simulations in two dimensions, performed with the code PLUTO and using the 
adaptive mesh refinement technique, are used to model interacting stellar and pulsar winds on scales ~ 80 times the 
distance between the stars. The hydrodynamical results suggest the suitable locations of sites for particle acceleration 
and non-thermal emission. 

Results. In addition to the shock formed towards the star, the shocked and unshocked components of the pulsar wind 
flowing away from the star terminate by means of additional strong shocks produced by the orbital motion. Strong 
instabilities lead to the development of turbulence and an effective two-wind mixing in both the leading and trailing 
sides of the interaction structure, which starts to merge with itself after one orbit. The adopted moderate pulsar-wind 
Lorentz factor already provides a good qualitative description of the phenomena involved in high-mass binaries with 
pulsars, and can capture important physical effects that would not appear in non-relativistic treatments. 
Conclusions. Simulations show that shocks, instabilities, and mass-loading yield efficient mass, momentum, and energy 
exchanges between the pulsar and the stellar winds. This renders a rapid increase in the entropy of the shocked structure, 
which will likely be disrupted on scales beyond the simulated ones. Several sites of particle acceleration and low- and 
high-energy emission can be identified. Doppler boosting will have significant and complex effects on radiation. 
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1. Introduction 

The binary sy stem PSR B1259-63/L S2883, consisting of 
a late O star (|Negueruela et al. 2011 ) and a 47 ms pulsar 
( Johnston et al.||1992p , is a powerful GeV and TeV emitter 
e.g. [Aharonian et al 



2005 



Abdo et al. 2011 Tam et al 



201l|. The physical processes behind the gamma-ray emis 



simTdetected in PSR B1259-63/LS2883 are believed to 
take place in the region where the stellar and pulsar winds 
collide, which would also produce X-rays and, once the in- 
teracting flows have left the system, radio emission (e.g., 
Tavani fc Arons||19971 [Johnston et aTlpMl [Dubusl|20M| 



Kirk||2000[ [Khangulyan et al.||2007[ |2011[ |2012[ ). The de- 
tails oFtEe dynamics of the interacting tiows are still poorly 
known, which affects the interpretation of the observed vari- 
ability and spectra not only in PSR B1259-63/LS2883 
but also in other gamma-ray emitting binaries that may 
host a non-accreting pu lsar (LS 5039 and LS I -1-61 303, 
see d i scussions in, e.g., [Chernyakova et al.||2006[ [Dubus 



l2006[ [Romero et al.||2007[ [Bosch-Ramon fc Khangulyaii 
2009 [Torres[[2011[[Bednarek[[2011[[Barkov fc Khangul^ 



201^ HESS J0632+057, see, e.g . [Bongiorno et al.j 



Moldon et al. 2011a). The processes underlying the non- 
thermal radiation in PSR B1259-63/LS2883 seem to be 
shocks that take place in the two-wind interaction region 
and accelerate electrons. These electrons most likely radi- 
ate from radio-to- X-rays through synchrotron and at GeV 
and TeV energies through inverse Compton emission (IC), 
although some gamma ray s may come from IC in the un- 
shocked pulsar wind (e.g., Tavani & Arons 1997 [Ball fc 



IFGL J1018.6— 5856, s ee, e.g., [Ackermann et al.J 
Abramowski et al.[|2012[ ). 



2011 



2012 



Hydrodynamical and magnetohydrodynamical, rela- 
tivistic and non-relativistic simulations of stellar/pulsar- 
wind int eractions have been c arrie d out in the past few 
years. In Romero et al. (2007) and Okazaki et al. (2011), 
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three-dimensional (3U), non-relativistic hydrodynamical 
(SPH) simulations were carried out to study, on binary 
scales but including orbital motion, the stellar/pulsar-wind 
int eraction in LS I +61 303 and PSR B1259-63/LS2883. 
In Okazaki et al. (2011), the authors explicitly studied 



the interaction with the equatorial disc present in the lat- 
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ter. Bogovalov et al. (2008 2012) performed axisymmctric, 
relativistic, hydrodynamical, and magnetohydrodynamical 
simulations of the stellar/pulsar-wind interaction, with the 
symmetry axis joining the two stars. The orbital motion was 
not accounted for, although the impact of the pulsar wind 
anisotropy was studied. Therefore, from previous work the 
flow evolution seems to be well-understood on scales up to a 
few pulsar/star separation distances (a) . One can conclude 
that the shape of the interaction region depends on the 
pulsar-to-star wind-momentum rate ratio (rj), the shocked 
pulsar wind reaccelerates when flowing away from the shock 
facing the star, and the bending of the interaction structure 
due to orbital motion is expected. 

An analytical study of the shocked flow evolution on 
scales larger than a indicated that Coriolis forces related to 
the orbital motion could amplify the shocked flow bend- 
ing, terminating the pulsar wind flowing away from the 
star with a strong shock, and enhance instabilit ies and 
mixing in the flow cont act discontinuity (CD) (Bosch- 
Ramon & Barkov 2011). The basic idea is that under 
typical ry- values, ~ 0.03 — 0.3, the kinetic luminosity of 
the pulsar wind dominates that of the star by a factor 
X ~ 30(?7/0.1)('(;w/2xl0^cm s~^)~^, where Ww is the stellar 
wind velocity. Such a high x-value, the interaction geom- 
etry, and the impact of the Coriolis force will eventually 
lead to the formation of a partially confined, high entropy 
flow with a strong pressure gradient outwards. This will 
lead to isotropization of the flow mass, momentum, and en- 
ergy fluxes, and to reexpansion, with the consequent loss of 
structure. In this context, non-relativistic two-dim ensional 
(2D) simulations done by Lamberts et al. (2012) showed 
that the shocked flow may keep a spiral-like shape up to 
large distances, although the adopted light-to-dense wind 
luminosity ratios, x "£ 1.25, and the non-relativistic nature 
of the simulations may prevent one from generalizing this 
conclusion. 

To study in detail the pulsar/star wind interaction on 
scales ^ a, plus a pulsar wind velocity u — > c, we per- 
formed 2D relativistic hydrodynamical simulations includ- 
ing orbital motion, where the flow is homogeneous and per- 
pendicular to the simulated plane. We adopted x = 30 and 
77 = 0.3, the latter resulting in a geometry of the interac- 
tion region similar, on spatial scales ~ a, to the 3D case 
for 77 = 0.1. For complementarity, we also explored the case 
with ;!^ = 60 and rj = 0.6. We note that for 77 < 1 and dis- 
tances to the star r ^ a, the pulsar wind gets diluted faster 
than the stellar wind in 3D than 2D. This means that in 
3D, the shocked winds flowing away from the system are 
deflected by the Coriolis force at a shorter distance than 
obtained in the present calculations. Moreover, instabilities 
are typically less disruptive in 2D. All this shows that the 
flows simulated here are likely to be more stable than in 3D 
simulations, which should be carried out in the near future. 

The paper is organised as follows: in Sect.[2]we present 
simulations of the orbital evolution of the flow structure; in 
Sect. [3] the impact of the value of the pulsar- wind Lorentz 
factor is studied, and the consequences of adopting a rela- 
tivistic or a non-relativistic approach discussed; finally, in 
Sect. HI we summarize our work and present a discussion of 
the implications of our results for the non-thermal processes 
in high-mass binaries hosting a young pulsar. 



2. Numerical simulations 

2.1. Numerical set-up 

The sir nulations were implern ented in 2D with the PL UTO 
codaM (iMignone et al. 2007 1, the pi ece-parabolic method 
(PPM) (ColcUa fc Woodward] pMf, an HLLC Riemann 
Solver (Mignonc & Bodo 2005), and using AMBrl through 
the Chombo codtf] ( .Colella et al. 2009^ . PLUTO is a mod- 
ular Godunov-type code entirely written in C intended 
mainly for astrophysical applications and high Mach num- 
ber flows in multiple spatial dimensions. The simulations 
were run through the MPl (message passing interface) li- 
brary in the cluster of Moscow State University Chebyshev. 

The simulated flows were approximated as an ideal, rel- 
ativistic adiabatic gas with no magnetic field, one particle 
species, and a polytropic index of 4/3. The adopted reso- 
lution was 96 X 96 cells and 6 levels of AMR, which gives 
an effective resolution of 6144 x 6144 cells. The size of the 
domain was x £ [0, 80 a] and y S [0, 80 a]. 

The stellar and pulsar winds were assumed to be 
isotropic and moving at a constant speed (see below). We 
thus assumed that the winds had already been accelerated, 
and did not consider the several effects affecting the state 
of the two winds before colliding: the ionization/heating ef- 
fects of X-rays from t he shocked wi nds on the stellar wind 
formation region (e.g. Blondin 19941; the pressure of the ra- 
diation from either the p ulsar or the shocked winds exerted 
on the stellar-wind (e.g. 



Gayley et al. 1997; 



of gr avity on slow and deiise material within the binary 



the effects 



(e.g. Blondin et al. 1991 



Okazaki et al. 



particularly relevant in Be and/or very c 



2011|), which are 
and 



lose systems 

the impact of the stellar radiation field on the pulsar-wind 
Lorentz fact or through C ompton braking (e.g. Ball & Kirk 
2000; Khangulyan et al.|[2007i ). We note that, even though 
all these factors may play a major role in particular sources 
and should be accounted for in more refined treatments of 
the problem, they can be considered at this stage higher 
order effects and have therefore been neglected. 

A pulsar-wind Lorentz factor F = 2 was adopted be- 
cause of resolution limitations. This i s smaller than the 

lO'' 



conventional value F 



2012 Aharonian et al. 2012 



- 10^ (see Khangulyan et al. 
and references therein), but 
high enough to reproduce the wind velocity contrast in the 
bow-like shock region, where instabilities trigger, and to 
capture important relativistic effects. Similarities and dif- 
ferences between the F-value adopted here, a more realistic 
one, and the non-relativistic case, are discussed in Sect. [31 
The orbit was assumed to be circular, with a period T — 
3 X 10^ s and radius a = 3 x 10^^ cm, and v^ = 0.01 c, which 
is typical for the isotropic/polar wind in massive stars. 
A sl ow equatorial wind is not considered here (see other- 
wise 



Okazaki et al.l 2011 



Bosch-Ramon & Barkov||2011|) 
„ values 



For the given 77, 1', and v^, values, the stellar and pulsar- 
wind momentum rates, normalized to the stellar mass-loss 
rate M-7 = (M/IO"'^ M© yr^^), are « 2 x lO^^'M.r and 
6 X 10^^M_7 g cm s~^, respectively, and the pulsar-wind 
luminosity « 9 x 10'^^ M_7 erg s~^. The numerical solution 
obtained can be used for other orbits with similar values of 



1 



Link 



http://plutocode.ph.unito.it/index.html 



^ Adaptive mesh refinement technique 
^ Link https://commons.lbl.gov/display/chombo/ 
'^ This effect is expected to be minor due to a star luminosity 
that is much higher than the spin-down one. 
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T /a. For high M- and small a-values, the shocked stellar 
wind may be radiative instead of adiabatic, enhan cing the 
growt h of the flow instability at the CD (see, e.g., [Pittard 



2009 



for two non-relativistic winds) 
At the beginning of the simulation, the unperturbed 
stellar and pulsar winds occupied the two grid halves, left 
and right, respectively. Both winds were supersonic, of neg- 
ligible temperature, and had a cylindrical symmetry around 
axes perpendicular to the simulation plane and crossing 
the wind sources. The given values of Ww and T, together 
with M and 77, determine the wind densities. The star 
and the pulsar were initially located at (39. 5o, 40o) and 
(40.5a, 40a), respectively, and the orbital motion, occurring 
in the plane XY , was set counter-clockwise. The simula- 
tions included neither the magnetic field nor an anisotropic 
pulsar wind, since their impact is expected to be small 



(Bogovalovet al. 2012|. 



2.2. Results 

Figures [T] and |4] show the evolution of the density map of 
the simulated flows during one orbital period, for rj — 0.3 
and 0.6, respectively. Physical times earlier than < 0.3 T are 
not shown since the structure had not yet reached a quasi- 
steady configuration on the scales of interest. The snapshots 
are obtained at t ~ 0.4 T, 0.7 T, and 0.9 T. For the last 
snapshot time, the tracer map is also shown in Figs. [2] (77 = 
0.3) and [5] (77 = 0.6), indicating how much pulsar (1) and 
stellar wind material (-1) are present in the computational 
cells. Finally, Figs.|3](r7 = 0.3) and|6](77 = 0.6) show a map 
of the module of thefour- velocity spatial component (three- 
velocity module times the flow Lorentz factor) . Arrows and 
their colours in all figures indicate the direction and module 
of the four- velocity spatial component, respectively. 

The flow behavior in the inner bow-shaped interaction 
regio n is very similar to that foun d in previous simulations 
(e.g., Bogovalov et al.|[2008 2012 1. In particular, the reac- 
celeration of the pulsar wmd after the termination shock 



Bogovalov et al. 2008) is also present in our calculations 
espite the adopted sraall F-value. This effect can be no- 
ticed in the velocity gradient around the shock apex shown 
in Figs. [3] and [6] with bulk Lorentz factors increasing from 
« 1.1 to 1.4 (see also Sect.lsl). On larger scales, however, 
the interaction structure starts to depart strongly from be- 
ing symmetric with respect to the two-star axis. Several 
orbital-motion effects are apparent. First, the stellar wind 
pushes clockwise the leading side of the interaction struc- 
ture, forming strong shocks in both the unterminated and 
the already terminated pulsar-wind zones. At the leading 
edge, and in the unterminated pulsar wind flowing away 
from the binary, a strong shock is formed at the distance of 
r ~ 5 a. On the trailing side, the termination takes place at 
around r ~ 10 a. Beyond these points, the shocked pulsar 
wind is deflected by an angle ^ 7r/2, and becomes strongly 
turbulent and increasingly loaded with stellar wind mate- 
rial. We note that, even for highly eccentric systems (such 
as PSR B1259-63/LS2883), the obtained solution stiU pro- 
vides a qualitative description of the shocked structure. The 
reason is that the characteristic flow speeds exceed the pul- 
sar orbital velocity by at least one order of magnitude. This 
implies, on the one hand, that the flow structure will adia- 
batically adjust to the changing a along the orbit. Moreover, 
the longer orbital timescale compared to that of the shocked 



flows will allow the latter to cover a significant fraction of 
a spiral turn before a has changed significantly. 

Turbulence is triggered by the Kelvin-Helmholtz (KH) 
instability in different regions of the two-wind shocked 
structure. In its leading edge, the instability grows at the 
two- wind CD, as also illustrated in Sect. [3j The instabil- 
ity couples with the perturbations generated at the pulsar- 
wind termination in the opposite direction to the star. At 
the trailing side, similar processes take place, although the 
instability grows substantially further, partially because the 
flow (re) terminates farther, partially because the flow lines 
diverge there. Bending will make the flow lines of the two 
shocked winds converge again, but it takes some time for 
the effect of the Coriolis force to propagate through the 
shocked structure and cause this convergence. 

The tracer maps in Figs. [2] and [5] indicate that already 
within the first orbital turn, i.e. on scales ~ (10 — 100) a, 
some amount of stellar wind has penetrated deeply into the 
shocked pulsar flow. This mass-load slows down the shocked 
pulsar wind, as shown in Fig. [Sj and in general makes the 
flow motion more chaotic, and rich both in weak shocks 
and turbulence. As seen in the figures, after one orbit the 
trailing side of the shocked stellar wind already approaches 
the leading edge, which has already become turbulent. 



3. On the pulsar-wind Lorentz factor 

As alluded in Sect. [Tl approximating the pulsar wind by 
a non-relativistic flow can underestimate the pulsar-to-star 
wind energy-flux ratio, for a given momentum-flux ratio, 
by a factor of v/2c. This can be (roughly) solved assum- 
ing a mildly relativistic pulsar-wind velocity while keeping 
the Newtonian treatment, which leads to moderately un- 
physical results. However, there are additional, more fun- 
damental, differences between relativistic and Newtonian 
flows that are relevant to the problem considered here. In 
this section, we discuss the impact of adopting a relativistic 
model for the simulated flows, and the influence of the spe- 
cific value of the pulsar-wind Lorentz factor. For instance, 
although the results obtained in the mildly and highly rela- 
tivistic cases are similar, the KH instability will tend to de- 
velop faster in the former. The impact of a non-relativistic 
approximation for the simulated flows is however stronger, 
since the shocked flow dynamics will be substantially dif- 
ferent when orbital motion is included, even for the mildly 
relativistic case adopted in Sect. [2j 

3.1. The colliding wind region at F = 10 

We simulated the central part of the interaction region as- 
suming a pulsar wind with F = 10 instead of 2. The ge- 
ometry and system properties were the same as in Sect. [2] 
including orbital motion, and rj = 0.3. The simulation was 
run long enough to reach a quasi-steady state on the con- 
sidered scales of both x e [0,4a] and y G [0,4a], with the 
optical star located at (2a, 2a). The adopted value of the 
wind Lorentz factor required an effective resolution of the 
computational grid of about three times the one adopted 
in Sect. [2j now with a domain size of 768 x 768 cells. The 
resulting maps for the distribution of density, temperature, 
and the spatial component of the four-velocity are shown 
in Figs. JTj [8J and [9J The plot shown in the bottom panel 
of Fig. fTTillustrates that the CD can change strongly on 
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Fig. 2. Tracer distribution by colour, where 1 is only pulsar 
wind and -1 only stellar wind, at i = 2.7 x 10^ s, for the 
case 77 — 0.3. 
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Fig. 3. Distribution by colour of the module of the spatial 
component of the four-velocity at < = 2.7 x 10^ s, for the 
case 1] = 0.3. 



timescales much shorter than the binary period. The re- 
sults of this numerical test show that strong shocks induced 
by the orbital motion are still present when the pulsar- 
wind Lorentz factor is higher, and the reacceleration of the 
shocked pulsar wind is stronger in the case of F = 10 than 
for r = 2. As shown in Sect. |3.3[ the shocks induced by 
the orbital motion are a general feature when adopting a 
relativistic flow model, and they should also be present in 
the ultra- relativistic case. Regarding the development of the 
KH instability at the CD, in the considered case the growth 
rate does not seem to be significantly reduced. However, 
this conclusion may not be generalized in a straightforward 
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way to the ultra-relativistic case, which is considered in the 
next section. 



3.2. Instability growth at the CD and F 

Although the shocked structure evolution is the result of 
the interplay of different mechanisms operating on differ- 
ent scales, the KH instability developing in the CD can by 
itself have a strong impact on the flow structure. To study 
the instability growth in the CD, we performed 2D test sim- 
ulations of the shocked pulsar and stellar winds adopting a 
mildly and an ultra-relativistic pulsar wind. The context of 
the simulations is the scenario described in Sect. [2] The cal- 
culations were performed with a relativistic hydrodynam- 
ical code whose characteristics are described for instance 
at the end of Sect. 2 in 



(|2C 
.3" 



2005). The size of 



Perucho et al 
the computational box was lO^"* cm x 3 
resolution, 1920 x 640 cells. The two shocked winds are ho- 
mogeneous in the direction perpendicular to the simulated 
plane, and are initially in thermal pressure balance. It was 
assumed that the shocked pulsar wind had already been 
reaccelerated, thus the thermal pressure was significantly 
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Fig. 7. Top: Density distribution by colour of the interac- 
tion region for a pulsar wind with F = 10 and r] = 0.3. The 
domain shown is limited to highlight the interaction region. 
Bottom: The same as in the top but few hours later, to illus- 
trate that the CD can change strongly on timescales much 
shorter than the binary period. 



lower than the ram pressure. The shocked stellar wind was 
located at the upper half of the computational box (see 
below) and was assumed to be initially at rest. 

The results for the density distribution illustrating the 
growth of seed perturbations are presented in Figs. 



10 and 



[TTJ The pulsar wind in Fig. [TO] has a bulk Lorentz factor of 
2, and a thermal pressure of about 10% of the ram pressure. 
In Figure 11 the pulsar-wind bulk Lorentz factor is 10, and 
the microscopic Lorentz factor in the flow frame, related to 
temperature, is 10^. This case is very different from the case 
considered in Sect. |3. 11 and would roughly correspond to a 




Fig. 9. The same as in Fig. [t] but for the module of the 
spatial component of the four-velocity. 



10 



and 



11 



shows 
of the 



shocked) F ~ 10^ . A comparison of Figs 
that, in the ultra-relativistic case, the development 
KH instability is significantly weaker, which is consistent 
with earlier numerical studies. In particular, a discussion 
of the reduction in the instability growth and related non- 
linear effects (namely mixing, disruption, an d deceleration) 
for high Lorentz-factor flows ca n be found in|Perucho et al. 



shocked and reaccelerated ultra-relativistic wind with (un- 



^^___ [Perucho et al. ( |2004b[ ), and [Perucho et al.| ( |2005f ^ 

in the context of planar and axisymmetric cylindrical jets. 
Although the situation considered in this study has been 
very much simplified, one may still draw the conclusion that 
the shocked structure in high-mass binaries hosting a pulsar 
is less likely to develop the KH instability than massive star 
binaries. However, even in the case of an ultra-relativistic 
flow, the figures show that there is still a significant growth 
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Fig. 10. Density distribution by colour in arbitrary units of 
the development of instabilities at the contact discontinuity 
after 2.8 (top) and 43 s (bottom) from the beginning of the 
simulation, for a mildly relativistic pulsar wind. The axes 
units are 1.7 x 10^" cm. 
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region within the system, where the orbital motion is neg- 
ligible, on scales > a the Coriolis force starts to deflect the 
flow, and inertia then becomes important. 

If one of the shocked winds is at least mildly relativistic, 
for it to be significantly deflected the exerted lateral pres- 
sure is to be a substantial fraction of the ram pressure of 
the flow. Since the plasma flowing out from the system is al- 
ready trans-sonic, this lateral pressure will be (F^x) higher 
than the thermal pressure, and the flow deflection will pro- 
duce shocks. These shocks, which are stronger for higher F, 
will trigger large CD perturbations, and can be a source of 
energy dissipation in the form of particle acceleration. 

In a Newtonian flow with the momentum and energy 
fluxes of the shocked pulsar wind, density and therefore 
inertia would be very small even if the plasma were rel- 
atively hot. Such a fluid could be smoothly deflected by 
sound waves, rendering a much more stable structure with- 
out much non-thermal activity. This conclusion also applies 
if the wind has a non-relativistic velocity but the same mo- 
mentum flux, and thus smaller energy flux, than its rela- 
tivistic counterpart. 

The effects discussed here balance the higher stability 
of the CD with an ultra-relativistic pulsar wind discussed 
in Sect. 
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Fig. 11. The same as in Fig. 



10 



but after 2.8 (top) and 34 s 
bulk Lorentz factor of the 



(bottom), for the case with a 
shocked pulsar wind of 10 (in the laboratory frame), and a 
microscopic Lorentz factor for the particles of 10^ (in the 
flow frame). 



of the KH instability on scales comparable to those of the 
binary system. This, plus the strong asymmetry of the in- 
teraction introduced by orbital motion, makes the genera- 
tion of shocks and turbulent flows at the CD very likely on 
scales of ~ a, regardless of the value of F. 

3.3. Relativistic vs Newtonian inertia 

Orbital motion introduces a strong non-stationary compo- 
nent in the dynamics of the shocked structure. In such a 
case, the inertia of the flow has a very strong influence, 
and there is a major difference between the relativistic and 
Newtonian case: the relativistic inertia is proportional to 
pressure when this is not negligible compared to the rest 
mass energy (as expected here), whereas Newtonian inertia 
is proportional to density. Therefore, although Newtonian 
simulations may reproduce the geometry of the interaction 
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60 



By itself, the instability growth rate would 
depend on (small) seed perturbations of the CD of some 
sort. However, under orbital motion, strong non-linear per- 
turbations are unavoidable, as shown in Sect. 3.1 We note 
that this conclusion is based on the relativistic nature of 
the flow, and does not apply to massive star binaries with 
colliding winds. 

4. Discussion 

The simulation results show that the colliding pulsar and 
stellar winds produce a coherent flow with low entropy up 
to only a few a. Beyond that point, for stellar and pulsar 
winds typical for gamma-ray binaries (possibly) hosting a 
non-accreting pulsar, the shocked flow becomes bent, non- 
ballistic, and highly turbulent, with strong wind mixing. 
It seems therefore unlikely that the shocked structure will 
remain as a spiral on scales ^ a (even for rj = 0.6 < 1). 
As seen by comparing the cases with rj — 0.3 and 0.6, the 
smaller the two-wind momentum flux ratio, the earlier the 
shocked structure is closed and disrupted by the orbital 
motion. Overall, the flo w evolution is quite similar to e arlier 
analytical predictions ( Bosch- Ramon fc Barkov||2011[ ). 

The adopted F-value already captures the main physics 
of the colliding wind region and the shocked structure 
formed there. Higher wind Lorentz factors may yield more 
stable flows at the CD in an idealized situation, but the 
impact of orbital motion enhances instability and triggers 
shocks, and overall the global picture should not change 
dramatically between F = 2 and F 3> 1. Newtonian simula- 
tions thus seem suitable only for the colliding wind region 
on its smallest scales, within the binary, because further 
out the impact of orbital motion will be much smoother 
than in the relativistic case, missing much of the insta- 
bility and shock generation. In addition, a non-relativistic 
pulsar wind, unlike its relativistic counterpart, could not 
re-accelerate and its speed could not then again approach 
c after being shocked. 

In addition to the pulsar-wind shock facing the star, 
and the unshocked pulsar wind itself, the present calcu- 
lations reveal additional sites of non-thermal activity (as 
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hinted already in Bosch- Ramon & Barkov 20111. Most of 
the pulsar-wind kinetic energy that is not reprocessed on 
the binary system scales is converted into internal energy 
farther out, either through strong shocks in the direction 
away from the star, or through weaker shocks and tur- 
bulence further downstream. Given their proximity to the 
binary system, the strong shocks behind the pulsar seem 
ideal sites for particle acceleration and high-energy emis- 
sion. For instance, farther from the binary, the expected 
weaker synchrotron and inverse Compton losses would al- 
low the accelerated particles to reach higher energies. This 
could explain the very energetic photons seen from LS 5039 
( Khangulyan et al. [2008a) . In addition, an emitter farther 
from the star woula allow gamma rays to avoid severe ab- 
sorp tion in the stellar photon field (e.g., Bosch- Ramon et al 
[2008 ) . Finally, the presence of emitting regions with differ- 
ent physical conditions may explain the broad-band spec 
trum in some gamma-ray binaries (e 
Zabalza et al., in preparation). 
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T^he complex hydrodynamics of the shocked flow on 
scales ^ a, with weak shocks and turbulence, indicates 
that radio-emitting electrons may be accelerated there in- 
stead of at the termination shock facing the star, where 
free-free radio absorption can be sever e, and the high en 



ergies of the accelerated electrons (e.g.. Kennel & Coroniti 
1984 1 are unsuitable for efficient radio emission. T'he radio 
electrons will initially move at mildly relativistic speeds, 
at least in most of the emitting flow, but for r > 100 a 
further mixing will likely slow down the flow. The spatial 
scale of the simulated flow corresponds to ^ 10 — 100 mil- 
liarcseconds at ^ 2 kpc distance, so the production and 
evolution of the radio-emitting electrons, and the structure 
of the underlying flo w, can be studied with current VLBI 



instrumentation (see Dhawan et ar]|2006 Ribo et aL]|2008 



[Moldon et a L]|2011a|b ) 

The mildly relativistic speeds downstream of the strong 
shocks in the pulsar wind imply that Doppler boosting 
can be signiflcant t here (as proposed in Khangulyan et al. 
see also, e.g., Dubus et al.|2010). However, we must 



2008b| 

emphasize that the motion of the whole interaction struc- 
ture is non-ballistic, with rather complex velocity distribu- 
tions both in terms of module and direction. This should 
induce through Doppler boosting complicated variability 
patterns with different timescales in the non-thermal emis- 
sion. This variability would appear on top of that intrinsic 
to the non-thermal activity itself. 
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